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Abstract 



We consider the dynamics of Rydberg states of the hydrogen atom driven by 
a microwave field of elliptical polarization, with a possible additional static 
electric field. We concentrate on the effect of a resonant weak field - whose 
frequency is close to the Kepler frequency of the electron around the nucleus 
- which essentially produces no ionization of the atom, but completely mixes 
the various states inside an hydrogenic manifold of fixed principal quantum 
number. For sufficiently small fields, a perturbative approach (both in classi- 
cal and quantum mechanics) is relevant. For some configurations of the fields, 
the classical secular motion (i.e. evolution in time of the elliptical electronic 
trajectory) is shown to be predominantly chaotic. Changing the orientation 
of the static field with respect to the polarization of the microwave field allows 
us to investigate the effect of generalized time-reversal symmetry breaking on 
the statistical properties of energy levels. 



I. INTRODUCTION 



Sixteen years ago, Bohigas, Giannoni and Schmit, in a seminal paper H, formulated 
the conjecture that quantum systems which are chaotic in the classical limit, generically 
have statistical properties of energy levels described by Random Matrix Theory (RMT) . 
While RMT had proven useful in characterizing nuclear spectra [[3], the conjecture was quite 
surprising since it claimed an applicability of RMT statistics to deterministic systems with 
very few degrees of freedom. While some counterexamples exist (see e.g. [|J), the conjecture 
remains widely accepted and played a very stimulating role in quantum chaos studies (see 
for reviews [§,§J). 

Depending on the symmetries of a given strongly chaotic system, the statistical proper- 
ties of its quantum spectrum fall into one of the three classes known from RMT: orthogonal, 
unitary and symplectic. The orthogonal class is associated with Hamiltonians invariant 
with respect to some anti-unitary symmetry. Typically this is true for time-reversal invari- 
ant systems, but also for more complicated anti-unitary symmetry (generalized time-reversal 
symmetry) as the product of a time-reversal with some discrete reflection. Roughly speak- 
ing, this happens if a basis can be chosen where the matrix elements of the Hamiltonian 
are all real. The corresponding statistical ensemble of (real symmetric) random matrices is 
referred to as the Gaussian Orthogonal Ensemble (GOE). In the absence of any anti- unitary 
symmetry, the corresponding class of complex hermitian matrices is known as the Gaussian 
Unitary Ensemble (GUE). Finally, specific considerations apply to half- integer spin systems 
with (generalized) time-reversal symmetry. Indeed, the energy levels of these systems are 
systematically two-fold degenerate (this is the well-known Kramers degeneracy). In the pres- 
ence of a geometrical symmetry like azimuthal symmetry, the Kramers degeneracy is hidden 
and the GOE statistics apply. With no geometrical symmetries, the Gaussian Symplectic 
Ensemble (GSE) has to be used. All these three classes of systems are characterized by level 
repulsion, that is zero probability of observing accidentally degenerate energy levels. This 
is to be contrasted with a generic behaviour of multidimensional integrable systems where 
close lying levels are uncorrelated and obey a Poisson statistics 0. For a generic Hamilto- 
nian system, where typically chaotic and regular motions coexist, the statistical properties of 
energy levels are not universal and are expected to be intermediate between the two limiting 
cases, Poisson and RMT. 

While the conjecture has been tested on a number of theoretical models (we refer the 
reader to reviews rather than numerous original papers), experimental verifications are 
much less abundant and restricted fully, as far as we know, to the orthogonal universality 
class ||. This is probably due to the fact that most experimental results are obtained 
for atomic and molecular spectra. In the presence of a uniform magnetic field, the time- 
reversal symmetry is broken, but a anti-unitary symmetry persists || . Thus, observing GUE 
statistics should require well controlled inhomogeneous fields on the atomic scale, which has 
not been achieved. 

An alternative is to study the eigenmodes not of the Schrodinger equation, but of a dif- 
ferent although similar wave equation. The best examples are two-dimensional microwave 
billiards (in fact three-dimensional flat cavities below the cut-off frequency) where the clas- 
sical Helmholtz wave equation is in fact equivalent to the time-independent Schrodinger 
equation. There, breaking of the time reversal invariance is possible by using magnetic de- 
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vices and GUE-type statistics have been observed experimentally [pUlQf. Still it is desirable 
to have at hand a quantum system where manifestations of the generalized time reversal 
symmetry breaking is experimentally accessible. The aim of this paper is to discuss an 
example of such a system. As mentioned above, breaking all anti-unitary symmetries using 
static fields is not possible. Hence, the idea is to use a time-dependent field acting on an 
atom. With an oscillatory field alone, the product of time reversal by symmetry with respect 
to one of the polarization axis is an anti-unitary symmetry. Hence, GUE statistics requires 
to combine a microwave field and a static field. 

The model we use, the Rydberg states of the hydrogen atom driven by a microwave 
field, has a long history of its own. The system attracted attention when it turned out 
that experimental results [[ill on the ionization probability as a function of the microwave 
field and frequency could be explained in terms of the underlying classical chaotic dynamics 
which results in a diffusive energy gain of the Rydberg electron and eventually to ionization 
||12|| . Many interesting phenomena (e.g. quantum localization) have been studied on this 
model both theoretically and experimentally (see e.g. review papers [p~3H T6|1 ) . The first 
experiments involved linearly polarized microwaves allowing for a rough description of the 
ionization thresholds via a simple one-dimensional time-dependent model where the motion 
of the electron is restricted to the microwave field axis |T7|. For other polarizations (available 



in recent experiments [18,19[), such a simple model is not possible: the simplest dynamics 



may be two-dimensional where the electronic motion is restricted to the polarization plane. 
Such models have been investigated both for circular polarization (CP) p0|-|22|| and general 
elliptical polarization (EP) |23|-p5 . 



When not only the ionization thresholds, but also the more subtle details of the dynamics, 
have to be studied, a full three dimensional atom must be considered. For LP microwave, 
the conservation of the angular momentum projection onto the polarization axis, L z , makes 
the dynamics effectively two-dimensional and time dependent. For CP case, while L z is not 
conserved, the transformation to the frame rotating with the microwave frequency removes 
the explicit oscillatory time-dependence, leading to a three-dimensional time independent 
problem. Both these simplifications are no longer possible in the general EP microwave 
field and the problem becomes truly three dimensional and time dependent, providing new 
challenges to the theory. It is this situation which we shall consider later on. 

However, instead of discussing typical "large" microwave amplitude for fields which lead 
to an efficient excitation of Rydberg atoms to higher excited states and finally to ioniza- 
tion, we shall consider a "weak" perturbation of an atom where the exchange of energy 
between the field and the atom is negligible. The fact that the Coulomb problem is highly 
degenerate provides us with the possibility, as shown below, to reach chaotic dynamics even 
for infinitesimally weak external fields. This is not in contradiction with the Kolmogorov- 
Arnold-Moser (KAM) theorem ]26| since the latter precisely does not apply for degenerate 
situations. Those weak fields do not really excite the atom, but rather couple and mix 
among each other the Uq states of a degenerate hydrogenic manifold with a given principal 
quantum number n . The degeneracy is lifted, states undergo small (but measurable) shifts 
that reflect the dynamics of the no manifold. Such a situation is referred to as intramani- 
fold dynamics (by contrast to a typical intermanifold coupling leading to atomic excitation 
and ionization). An intramanifold chaotic behaviour is quite attractive from the theoretical 
point of view. It yields an effective quantum Hamiltonian acting in the n space, whose 
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eigenvalues are the energy shifts, and is represented by a finite dimensional matrix of size 
n 2 , (modulo the remaining point symmetries) with no cutoff errors. The semiclassical limit 
is realized by letting n — > oo. 

The first attempt to produce some intramanifold chaos used the hydrogen atom in uni- 



form crossed magnetic and electric fields pT|,|2"8||, which was later extended to arbitrary 



mutual orientations of the two fields p9 |. The authors, exploring second order perturbation 
theory, have observed signatures of intramanifold chaos in the quantum spectrum. However, 
the situation is somewhat complex: the first order term (in the two external fields) is always 
integrable; it is only when combined with a second order term of a comparable magnitude 
that some noticeable region of chaos can be created. But this implies the application of 
large fields, and higher order terms in the perturbative expansion are of importance. That, 
unfortunately, has been neglected by the authors. Moreover, the high electric field values 
used lead then to an extremely fast field ionization which blurs any long-time effects in the 
dynamics. Especially, there are no more bound states in these conditions, but only broad 
overlapping resonances, which implies that the statistical properties cannot be described by 
one of the three standard ensembles of random matrices |30[ . 



Recently, we have proposed |Jl]] another system revealing chaotic intramanifold dynamics 
which does not suffer from the mentioned deficiency. This is an hydrogen atom driven 
resonantly (i.e. with the frequency u which is an integer multiple of the Kepler frequency 
cuk = l/ n o) ^ an elliptically polarized microwave with the possible addition of some static 
weak electric field. The latter allows us to break any generalized time-reversal symmetry of 
the system. 

Since the driving field is periodic, by applying the Floquet theorem, one can find the 
eigenstates of the system (the so-called Floquet or dressed states). The eigenenergies of the 
system are referred to as quasienergies of the system and are defined modulo Uu. As we have 



discussed shortly for a 2 : 1 resonance pl| , the statistical properties of the quasienergies 
reveal convincingly the symmetry breaking. Here we discuss the same system, i.e. an 
hydrogen atom driven by a resonant EP microwave field, both in the presence and in the 
absence of the static field. 

The paper is organized as follows. In Sec. II we derive the quantum perturbative Hamil- 
tonian for our system and its semiclassical counterpart. They are used in Sec. Ill to analyze 
the behavior of the atom in a pure microwave field and in Sec. IV for a microwave field 
combined with a static electric field. The summary and the future perspectives form the 
content of the concluding section. 



II. PERTURBATION APPROACH 

We consider a realistic three-dimensional hydrogen atom placed in a static electric field 
and driven by an elliptically polarized microwave field. We define the z-axis as perpendicular 
to the plane of the polarization of the microwave field. The Hamiltonian of the system, 
in atomic units, neglecting relativistic effects, assuming infinite mass of the nucleus, and 
employing dipole approximation reads: 

P 2 1 

H=^- h F(xcosuut + ay smut) + E • r, (2.1) 
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where F, a and 00 stand, respectively, for the amplitude, degree of elliptical polarization and 
frequency of the microwave field while E denotes the static electric field. 

As mentioned in the Introduction, we are interested in effects due to weak external 
fields. In the absence of any external field, the energy levels of the system are — 1/2wq with 
no the principal quantum number; the degeneracy of the no hydrogenic manifold is Uq. Even 
very weak fields mix strongly states within the manifold; on the other hand, it is perfectly 
justified to treat coupling to other manifolds perturbatively. We shall do so first quantum- 
mechanically with the help of the effective Hamiltonian approach ||32|| . Then, we construct 
its (semi-)classical equivalent. The latter allows us to study the character of the classical 
motion and search for parameters corresponding to chaotic dynamics. 



A. Quantum perturbation method 

For any time-periodic Hamiltonian, the Floquet theorem states that the most general 
solution of the Schrodinger equation can be written as a linear combination of the "Floquet 
states", which are time-periodic functions and eigenstates of the so-called Floquet Hamilto- 



nian 7iFioquet @ of the system: 



H Floquet I <£(*)> = (h - i^J I 0(t)> = e I c/>(t)) (2.2) 

where e and | </>(t)) are respectively the quasienergy and the Floquet state. The solu- 
tions of Eq. ( ^.2[ ) have to satisfy the usual boundary conditions in configuration space and 
periodic boundary condition in the time coordinate (by construction, Floquet states are 
time-periodic). Hence, the Floquet Hamiltonian 7^Fioquet acts on an extended Hilbert space 
containing also the time coordinate. 

For studying the quasi-energy spectrum, is it convenient to choose, as a basis of the 
atomic Hilbert space, the Sturmian functions (see e.g. p4|-|36| for details of this application) 
I n, L, M)( A \ where L and M are the usual angular and magnetic quantum numbers respec- 
tively while n > L + 1 labels the radial functions whose number of nodes is n — L — 1. A is 
a scaling parameter (unit of length in configuration space) for the Sturmian functions. For 
A = no, the Sturmian functions with n = no are the exact hydrogenic states, eigenstates of 
the unperturbed atom. As we intend to describe the dynamics inside the no-manifold of the 
hydrogen atom, we choose to keep A = n in all calculations. 

Along the time coordinate, we choose the usual oscillating exponential functions as a 
basis of time-periodic functions. They are labeled with an integer index K and defined by: 

(t\K) = e~ iKuJt (2.3) 

The whole Hilbert space of the atom + periodic perturbation system (Floquet Hamilto- 
nian) is spanned by the tensor product of the configuration space and time basis: 

I n, L, M,K)=\ n, L, M) ( " o) ® \K). (2.4) 

In the dressed atom language, K may be loosely identified with the number of photons 
exchanged between the atom and the field. 
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The Sturmian basis is not orthogonal, but it satisfies the following relation: 



£ \n,L,M,K)(n,L,M,K\±- = l (2.5) 



n,L,M,K 



The advantage is that, when written in the Sturmian basis, all the matrix elements repre- 
senting the various parts of the Floquet Hamiltonian have some strong selection rules. The 
selection rules on K trivially come from the Fourier expansion of the time dependences; the 
selection rules on L and M originate from the angular dependence of the various operators. 
The selection rules on n are far from obvious and are at the heart of the definition and 
properties of the Sturmian functions [p7f| . 



If we define: 

U = F(xcosut + ay sin oot), 

V = E-r, (2.6) 
we obtain the following selection rules for the matrices representing these operators in the 



Sturmian basis : 



An = 0, ±1, AL = 0, AM = 0, AK = 0, for H and unity operators, 

An = 0,±l,±2 AL = ±1, AM = 0, ±1 AK = 0, fori/, (2.7) 

An = 0, ±1, ±2, AL = ±1, AM = ±1, AK = ±1, for U. 

In addition, all matrix elements are known in closed forms and involve only square 
roots of rational functions of the various quantum numbers. Note that, because of the non 
orthogonal character of the Sturmian basis, calculating the Floquet quasi-energies requires 
to solve a generalized eigenvalue problem rather than a standard one. This is the price to 
pay for getting sparse matrices. 

The exact calculation of the quasi-energies is possible only numerically. However, we 
are interested in the situation where both the static and the microwave fields are weak. 
Thus, a perturbative expansion is convenient. We will now perform it at the lowest non- 
vanishing order for each external field. Because, the zeroth order eigenstates are highly 
degenerate, we have to use degenerate perturbation theory. A convenient formulation is to 
use an effective Hamiltonian which, at any order of the calculation, has the same spectrum 
as the initial Hamiltonian, but acts only inside the degenerate hydrogenic manifold. The 
details of the method are given in |32| . At first order, the calculation is trivial and the 



effective Hamiltonian is just the projection of the perturbation onto the manifold we 
are considering. If P denotes the projector onto the degenerate (n , K = 0) manifold, it is 
simply: 

H (1) = P(U + V)P = PVP (2.8) 

since U always changes K by one unit. The non- zero matrix elements of are those of 
V with An = 0. Thus is proportional to the static electric field. 

The lowest non-vanishing contribution of the microwave field is at second order. It has 
the following well-known formal expression: 
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= PUQ- ttQUP (2.9) 

where A = — l/2 n o * s the unperturbed energy of the hydrogenic manifold and Q = 1 — P 
is the projector onto the subspace complementary to the hydrogenic manifold. 

Explicit calculation of is not straightforward. Indeed, if one expands the 1/(E — H ) 
onto the eigenstates of Ho, one obtains a infinite sum over the discrete states and continuum 
of the atomic spectrum. The trick is to use the non-orthogonal Sturmian basis defined above. 
Indeed, the projector P has the following simple expression is this basis: 

P=E \n ,L,M,0){n ,L,M,0\^-, (2.10) 

L,M Zr 

and consequently 

Q= \n,L,M,K){n,L,M,K\± (2.11) 

L,M,(n,K)^(n ,0) 

The last step is to calculate the matrix element of the l/(-Eo — Ho) operator in the 
Sturmian basis. It is simply accomplished by noting that the operator (Eq — Hq) is diagonal 
in L,M and K and tridiagonal in n (i.e. connects only state n to states n — l,n and n + 1). 
Thus the matrix elements of 1/(E — H ) are simply obtained by solving a triadiagonal set 
of coupled equations in each (L, M, K) subspace coupled to the initial state. 

Finally the whole effective quantum Hamiltonian inside the no-manifold reads 

H cS = H {1) +H {2 \ (2.12) 

This Hamiltonian takes into account the direct coupling between the levels due to the pres- 
ence of the static electric field (the term proportional to E) and the indirect coupling through 
all levels of other manifolds, i.e. process of absorption and emission of microwave photons 
(the term proportional to F 2 ) [|S2]] . Because of the selection rules on the U and V operators 
and the simple algebraic structure of the effective Hamiltonian, H c g itself has the following 
selections rules: An = 0, AK = (by construction) and AL = 0, ±1, ±2, AM = 0, ±1, ±2. 
The diagonalization of H c g, i.e. of the sparse banded matrix of dimension n^, by standard 
routines, yields quasienergies of the system. 

The method has been tested in limiting situations, e.g., for parallel weak static and 
linearly polarized fields where quasienergies resulting from the effective Hamiltonian could 
be compared with exact diagonalization values. The results for EP presented below are 
obtained for field amplitudes for which excellent agreement between exact and perturbative 
results have been found in the limiting cases. 



B. Semiclassical perturbation method 

The general prescription for the semiclassical quantization of a time-periodic system has 
been described in |39|| . Recently, we used a similar procedure for a LP microwave [M, for 
a static electric field parallel to a LP microwave [[lip and for a two-dimensional model atom 
in the EP microwave case [p4] , p5[| . 
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The method requires first to define a classical Hamilton function for which one can use 
the usual semiclassical quantization rules. This is done by passing to the extended phase 
space, defining the momentum p t conjugate to the t (time) variable. It yields the new 
classical Hamiltonian, 7^Fioquet = H + p t f26|, being the classical analog of the quantum 
Floquet Hamiltonian, Eq. fl2.2|). 

As the next step, we express the Hamiltonian in action-angle variables of the unperturbed 
Coulomb problem Due to its high symmetry, several choices are possible. The 

standard solution is to consider the canonically conjugate pairs (J, 0), (L, \&) and (M, $). 
J is the principal action (corresponding to the principal quantum number), i.e. the total 
action along an unperturbed Kepler elliptical trajectory of the electron. It is simply related 
to the size of the ellipse. The corresponding angle, G, determines the position of the electron 
on its elliptical trajectory and depends linearly on time, G ~ u>Kt, for an unperturbed atom. 
L is the angular momentum, \l/ the angle of rotation around the axis defined by the angular 
momentum vector. Similarly, M, $ denote the projection of the angular momentum on 
the laboratory z-axis and the angle of rotation around that axis, respectively. The shape 
of the ellipse is best described by its eccentricity e - 



'1 — L 2 / J 2 while its orientation in 
the configuration space is determined by the Euler angles ($, arccos M/L, \l/) as defined by 



Goldstein 44 



Using these canonical coordinates, it is possible to write down the full Floquet Hamilto- 
nian. We now specialize to the resonant case where the microwave frequency is an integer 
multiple of the Kepler frequency of the unperturbed electron, i.e. u = uj/ujk = m. The 
corresponding action is: 



J = n 



CO 



-1/3 



A' 



m, 



(2.13) 



n is interpreted as the principal quantum number of the quantum hydrogenic manifold 
where the resonance takes place. 

In the absence of any external field, the variables (J, L, M, $) are all constant while 
Q evolves linearly in time (see above). Hence, in the presence of weak external fields, the 
motion along the variable will be much faster than along the other coordinates and the 
secular perturbation theory can be used: it averages over the nonresonant terms and 
yields the approximate resonant Hamiltonian of the form 



where 
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In 2 



3m 2 
2^ 



J 2 + FT m cos(6 - 5) + £7 + p t 



(2.14) 



--Tin 



M 

cos ip sin 6 ( cos $ cos ^ — — sin $ sin \1/ 

M 



+ sin if sin 9 ysin $ cos ^ + — cos $ sin \1/ 

/ M 2 
+ cos 6]ll — — sin \1> 




(2.15) 



and 
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= m6 - cut, J = LJ^. p t = Pt + —, (2.16) 

m m 

The secular variables O, J,p t are slowly varying variables obtained by substracting the un- 
perturbed resonant quantities. represents the phase drift of the electron along the elliptical 
orbit and J the distance (in action) to the exact resonance. 

The orientation of the static field with respect to the z-axis is determined by the usual 
spherical angles, ip and 9. The expression for 7 looks complicated, but it is actually 
nothing but the component of the average atomic dipole on the static field axis. Simi- 
larly, r m (L, M, $; a) and S(L, M, <£>; a) just represent the amplitude and the phase 
of the atomic dipole at the microwave frequency. They can be obtained simply from the 
Fourier components of the electron position on an unperturbed elliptical trajectory. The 
explicit, rather complicated formulae for 8(L, ty, M, a) and T m (L, M, $; a) are given 
by Eq. (2.15) and Eq. (2.16) of JEJ, respectively and are reproduced in Appendix A for the 
convenience of the reader. 

The last stage is to quantize the system using the approximate Hamiltonian, Eq. ( p. 14 ). 



As any explicit time dependence has disappeared in the secular Hamiltonian, the quantiza- 
tion of pt is trivial. Taking into account that eigenstates of the Floquet Hamiltonian have to 



be time-periodic, this yields p t = ktu (where k is an integer number) |39l , f40| which ensures 
the periodicity of the quasienergy spectrum with a period u>. 

The radial motion, i.e. in the (J, 0) effectively decouples from the secular motion of 
the elliptical trajectory, i.e. in the (L, \I/,M, <3>) space f4~2"| , |4"5| . While considering the ra- 
dial motion, the effective hamiltonian for the secular part is approximately constant (for 
a detailed discussion as well as possible counterexamples in some cases see [0). In effect 
the quantization resembles in spirit the Born-Oppenheimer approximation. One may first 
quantize the radial motion keeping the secular motion frozen. The radial motion exhibits a 
pendulum-like dynamics whose quantum eigenvalues are given by the solutions of the Math- 
ieu equation (see the similar treatment for one-dimensional systems |filf-|4"9[1). Hence, we can 



define an effective Hamiltonian acting in a reduced (L,ip } M, 0) phase space just replacing 
the (J, 0) part of H scc by the quantized energy levels of the pendulum. In this paper, we are 
interested in the ground state of the pendulum, thus, the quantization of the radial motion 
yields the following effective Hamiltonian for the secular motion: 

3m 2 

H eS = -— T a {q)+E 1 , (2.17) 

where the constant terms —l/2n^ and ku are omitted (7i e fr denotes the shift from the 
unperturbed energy level of the atom), and: 

q = f^r m . (2.18) 



is a dimensionless parameter. ao(q) is the Mathieu parameter |5(J corresponding to the 
ground state of the pendulum. We can introduce scaled quantities: 

F = n 4 F (2.19) 

E = n%E (2.20) 
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L = — 
n 

M 

M = — 

n 

r 

■p 771 

1 771,0 — ~ 

70 



L 
— 


(2.21) 


no 


M 


(2.22) 




n 


r 


(2.23) 






_7 


(2.24) 


ng' 





The effective Hamiltonian is: 



Ties = - JV«o(<?) + 5r7o- (2.25) 
8«o ™o 

For a large microwave field amplitude or in the deep semiclassical limit, i.e. n — > oo, 
we may employ the asymptotic expression, for large q, of the Mathieu parameter |5(J 



ao(q) -> -2g + 2^. (2.26) 

This corresponds to the case where the pendulum is localized near its stable equilibrium 
position (whose energy is —2q), its zero-point energy in the ground state being calculated 
in the harmonic approximation (hence the I^Jq term). 

In the opposite limit, i.e. for g C 1, another approximation exists ao(q) ~ — q 2 /2 [50 . 
This corresponds basically to a very weak trapping pendulum potential, where the motion 
is basically the free motion only slightly perturbed (at second order in the energy) by the 
potential. This is the classical counterpart of the quantum perturbative approach developed 
above, the equivalent of the "no n-mixing" approximation. In the following, we restrict 
ourselves to this case as more appropriate for moderate hq and very small Fq. Hence, the 
final expression for the effective secular Hamiltonian we are going to deal with is 

F 2 E 

HeS = ^2 F ™.o + — To- (2-27) 

This Hamiltonian is a semiclassical counterpart of the quantum effective Hamiltonian, 
Eq. Q2.12p , namely first order in the static electric field and second order in the resonant 
microwave field. 

To calculate the quasienergies semiclassically, one should quantize the secular motion 
determined by the Hamiltonian ( p. 27] ). It has been done in simpler situations (e.g. an 
hydrogen atom perturbed by a linearly polarized microwave field and a parallel static electric 
field [0). Then, the secular motion is integrable and its quantization straightforward 
using the WKB quantization rule. The present secular problem has two degrees of freedom 
and turns out to be non-integrable except for some limiting cases. Therefore, a detailed 
analysis of the classical motion in the phase space of secular variables is necessary for possible 
comparison with quantal data. 

III. PURE MICROWAVE PERTURBATION CASE 

Let us consider first the perturbation of an hydrogen atom by an elliptically polarized 
microwave field in the absence of the static field. Our previous studies of resonant driving of 
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the atom were restricted to the simplified two-dimensional model atom where the electronic 
motion is restricted to the polarization plane [p^ , Pq| . Then, the classical secular motion is 
one-dimensional and application of the WKB quantization rule gives very accurate results for 
quasienergies of the system. In a realistic three-dimensional model, a similar procedure is no 
longer possible as the effective classical Hamiltonian for a general EP field is not integrable. 
The secular motion of the atom for weak microwaves, see Eq. (|2.27 ), is determined by the 
Hamiltonian 

Kfr = 3^r^ j0 (L , *, M , $; a). (3.1) 

The integrable motion is obtained in the limiting polarization cases, i.e. a = or 1, only. 
That is, for the LP field, the angular momentum projection on the polarization axis is a 
constant of motion. For the CP case, because of circular symmetry of the field, $ becomes 
a cyclic variable and the secular motion is also effectively one-dimensional. Clearly, for a 
close to these limiting values, the secular motion will be close to integrable. With this in 
mind, looking for chaotic dynamics, we take a = 0.4 in the following (we have verified that 
this value is representative for a typical EP behaviour). 

Eq. ( p.l|) shows that the structure of the classical phase space of the secular motion 
depends only on the value of 7Y c fj/-F 2 (beside the integer m labeling the primary resonance). 
In other words, if the secular motion is non-integrable for some finite field amplitude, it re- 
mains non-integrable even for infinitesimally small amplitude. This clearly demonstrates the 
inapplicability of the Kolmogorov- Arnold- Moser theorem to the highly degenerate Coulomb 
problem. On the other hand, the time scale of precession of an electronic ellipse is affected 
by the strength of the perturbation; for very small F , it will be extremely slow, but the 
trajectories of the secular motion do not depend on F . 

Consider the principal 1:1 resonance case, i.e. m = 1. To focus on the phase space 
structure, we have plotted Poincare surfaces of section (SOS) for a few values of H e g/ Fq in 
Fig. |l|. For high values of the secular energy, the motion is generally regular. However, for 
an energy interval around H-ck/Fq = 0.3, the mixed character of the motion is apparent, a 
quite large chaotic layer is clearly visible. 

Switching to the 2:1 resonance case (right hand side of Fig. |I|) we find, as previously, 
regular phase space structures for high energy and mixed character of the motion for lower 
energies. By visual inspection, the secular motion for 2 : 1 resonance looks "more chaotic" 
with smaller regular islands embedded in a pronounced chaotic layer. Although we have 
searched quite extensively, we could not find values of a and the energy corresponding to 
fully chaotic motion. Always at best tiny regular islands have been found. 

The mixed character of the secular motion should have consequences on the statistical 
properties of the quasienergy spectrum inside the hq manifold. As the system possesses anti- 
unitary symmetry invariance, i.e. is invariant under time-reversal combined with the y — > — y 
transformation, see Eq. ( |2.1| ), the statistical properties are expected to reflect intermediate 
behavior between the Poisson and GOE character. To calculate the quasienergy spectrum, 
we employ the quantum effective Hamiltonian Eq. ( |2.12| ). One should take care of discrete 



symmetries of the system. That is, the system is invariant under the z — > —z transformation 
as well as the parity combined with translation in time by tt/uj transformations. Thus, the 
whole spectrum of the n -manifold splits into four independent smaller spectra which are 
unfolded independently in order to study level statistics. 
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The dynamics of the levels belonging to the no = 20 manifold as a function of the 
polarization degree a, are shown in Fig. [|, for the 1:1 and 2:1 resonance cases. In each panel, 
for clarity, there is only one of the four independent sub-spectra plotted. Qualitatively, the 
level dynamics reflects the character of classical motion, i.e. for high energies, one cannot 
see level repulsion and there are apparently level crossing (actually small avoided crossings). 
At lower energy, the level dynamics is more irregular with plenty of avoided crossings, a 
clear signature of classical chaos in the system, compare Fig. ^ which shows this region in 
more detail. 

To make the comparison more quantitative, we have calculated the cumulative nearest 
neighbor spacing (NNS) distributions, taking levels in the energy intervals H c s/Fq 6 [0.02 — 
0.045] and [0.00035 — 0.0018] for the 1:1 and 2:1 resonance cases, respectively. The intervals 
have been chosen to correspond to chaotic behavior as much as possible. Fig. ^ presents 
the results for the 1:1 resonance for principal quantum number n around 55 and for about 
twice bigger no, i.e. around 100. The similar results corresponding to the 2:1 resonance are 
plotted in Fig. [5[ In all cases, one can observe that the numerical data are intermediate 
between the Poisson and GOE statistics. However, the behavior of the 1:1 case is closer to 
Poisson while the 2:1 one to GOE character, in agreement with the more classically chaotic 
behaviour in the latter case. 

Quantitative measures can be obtained by fitting the data to some theoretical NNS 
distributions. Quite often, this is being done directly on the histograms for the NNS distri- 
butions. Such a procedure is, however, bin size dependent and should be avoided. We prefer 
to fit the integrated (cumulative) distributions which do not suffer from this ambiguity. 
There are several possible choices for the theoretical distributions. Berry-Robnik statistics 
51] corresponds to a superposition of independent Poisson and GOE spectra with relative 



weight q corresponding to the relative volumes of regular and chaotic parts in the classical 
phase space. Others possibilities are the phenomenological Brody and Izrailev ]53|] dis- 
tributions. The explicit expressions for the distributions used can be found in Appendix B. 
Although the Berry-Robnik distribution relies on some reasonable theoretical grounds while 
the two other ones are purely phenomenological, it is commonly accepted that, for not very 
small effective h, the Brody distribution - or much less known Izrailev one - works better 
than the Berry-Robnik proposition. The latter works well only in the very deep semiclassical 
limit (very small effective H) where tunneling between regular and chaotic parts of the phase 
space is negligible. For lower lying states, the "regular" and "irregular" part of the spectrum 
interact via tunneling, leading to level repulsion between states in the two groups. 

This behaviour is clearly visible in Figs. (| and |5|, (b) and (d). There, the Berry-Robnik 
distribution predicts much more small spacings than actually observed. On the other hand, 
except for small spacings we have found that Berry-Robnik statistics fits best in most of the 
cases, compare Fig. [| and Fig. [|. 

The obtained fitted values of the Berry-Robnik parameter and the parameters for the 
Brody and Izrailev distributions are given in Table |. While, for the latter cases, the param- 
eters have little physical meaning, as mentioned above, q in the Berry-Robnik distribution 
should measure the fraction of the chaotic phase space. Qualitatively, the value obtained 
agrees well with classical SOS plots. 

We have also studied the spectral rigidity, i.e. the A 3 statistics [@,|]||, in order to obtain 
some information on the long range correlations in the spectra of our system. Spectral 
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rigidity gives an independent information about the relative measure of the chaotic part of 
the classical phase space. For a superposition of independent Poisson and GOE spectra, one 
obtains |54| 



A 3 = A 3 P ™((1 - q)L) + Ai OE (qL). (3.2) 

We have fitted our numerical results with this distribution. The results are plotted in Fig. |B| 
while values of the fitted parameter are put in Table |. It is well known that, at large L, 
the spectral rigidity deviates from any universal behaviour and saturates. This is a non- 
semiclassical effect which should take place for larger and larger L as the effective h goes to 
zero. To fit the parameter q, we have taken into account only data up to L = 10 in order 
not to enter the saturation regions visible in Fig. |6|. 

Comparing the Poincare SOS, NNS distributions and A 3 statistics of the calculated data, 
the qualitative agreement between the classical dynamics and the corresponding quantum 
statistical properties is apparent. The fitted values for the relative measure of the chaotic 
part of the phase space coming either from the NNS or the A3 distributions agree perfectly 
and match well the visual aspect of the SOS. 

The NNS distributions change a little when the principal quantum number is modified. 
The n ph 55 case reveals slightly stronger level repulsion (and, therefore, a "more chaotic" 
character) than the no « 100 case. This suggests that some tiny regular structures in the 
classical phase space are not resolved for no = 55 but are for uq = 100. The long range 
correlations are more dramatically sensitive to n Q . The saturation of A 3 starts at about 
twice larger distance in L for the higher n value, see Fig. || This is in agreement with the 
theory of Berry ]55| where the critical L value scales as 1/K (the effective Planck constant 



in our problem is 1/no). For large L, the A3 statistics saturates at almost twice bigger 
value for no ~ 100 than no ~ 55 again in a qualitative agreement with theory. The latter 
[ 55H predicts (A 3 (oo) ~ 1/h) for a regular spectrum and (A 3 (oo) ~ ln(l/?i)) for a strongly 



chaotic system. We deal with an intermediate system with mixed phase space, thus we 
expect the numerical A 3 (oo) value to lie in between the two limits. This is indeed the case. 

Finally, let us briefly argue what happens for larger microwave field amplitudes. The 
secular motion, in the weak field limit, is determined by the effective Hamiltonian, Eq. ( |3.1| ). 
which in turn is a function of T m> o- Increasing the field amplitude one leaves the validity 
range of the Eq. ( |3.1|) and enters the region where Eq. Q2.26Q is applicable. Then the orbital 
electronic motion becomes localized inside a resonance island. The secular motion, however, 
remains unchanged because the new effective Hamiltonian is again a function of r m only. 
This means that the spectral statistical properties for higher field amplitude (of course not 
so big as to produce strong intermanifold mixing) will be the same as the ones presented 
here. For the high secular energy the motion will be also regular. Thus, one may expect 
that the nonspreading wavepackets predicted using the two-dimensional model p4| , p5| will 
also exist in the real three-dimensional world. 



IV. MICROWAVE PLUS STATIC ELECTRIC FIELD 

In this section, we discuss the intramanifold behaviour of the hydrogen atom exposed to 
a resonant microwave field and a static electric field of arbitrary mutual orientation. 
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For small field amplitudes, the effective secular Hamiltonian is given by Eq. ( |2.27| ). The 
Hamiltonian is the sum of two terms - the first one proportional to Fq (square of the scaled 
microwave field), the second one to Eq (scaled static field). For arbitrary mutual orientations 
of the two fields, the two terms have incompatible symmetry properties and, when having 
comparable magnitudes, induce a globally chaotic behavior. 

Eq. (|2.27|) has some well defined scaling properties with the field strengths F and E . 
Let us define the reduced microwave strength 

F=Fy /E = F 2 n 6 /E (4.1) 

and the reduced Hamiltonian 

H = H eS n 2 /E = -^fl,o + 7o- (4.2) 

The classical phase space structure depends only on the value of Ti and T (beside the static 
field vector orientation and the polarization of the microwave field), but not on the detailed 
values of no, Eq, F and the secular energy. Of course, weaker fields imply a slower secular 
motion, but this does not affect the structure of the phase space. In the quantum mechanical 
picture, the energy splitting of a degenerate hydrogenic manifold also depends on absolute 
values of the fields, but the structure of the levels does not. 

The application of the static electric field allows us to break any anti-unitary symmetry 
of the system. It is only for E x = (i.e. (p — 7r/2) or E y — (i.e. <p = 0) that the system 
has some anti-unitary invariance, under the combination of the time-reversal transformation 
with reflection with respect to the yz or xz plane, see Eq. (f2.1|). 

In a previous letter we have considered the system in the case of the 2:1 resonance 



driving, i.e. when the microwave frequency is twice the Kepler frequency. Let us consider 
here first the principal 1:1 resonance, i.e. m = 1. A possible signature of the anti- unitary 
symmetry breaking would be to observe level repulsion with stronger than linear repulsion 
(i.e. with P(s) oc s@ with (3 > 1 for small s). Clearly, it is desirable to have a predominantly 
chaotic classical dynamics, as a transition from GOE to GUE statistics is expected. To this 
end, we have to find values of the fields parameters which maximize chaoticity of the system. 
After rather extensive searches, we have found that a = 0.4, T = 10 and 9 = 7r/4 are a good 
choice for that purpose. The remaining spherical angle (p determining the orientation of the 
static electric field vector, is used to control the breaking of the anti-unitary symmetry. 

In Fig. [7|, we show Poincare SOS for a few values of Ti, for two cases: (p = and 7r/2. 
One can see that, for <p w tt/2, a predominantly chaotic structure space exists in phase space 
for a large range of secular energies. For ip = tt/2, the generalized time reversal invariance 
holds, as mentioned above. Thus, to study the symmetry breaking, it is interesting to e.g. 
decrease (or increase) (p gradually, collecting quantum data for some ip values. Note that 
the addition of the static field tends to make the classical dynamics slightly more chaotic, 
compare Figs [I] and [?| 

The data are then analyzed as in the previous Section - an example is shown in Fig. |8| 
for two different orientations of the electric field <p = {0.47T, 7r/2}. For each (p, data have 
been collected for principal quantum number no in the range 50-59. Then, we have chosen 
levels from the scaled energy interval H G (0, 1.4), unfolded each spectrum and calculated 
NNS distributions and spectral rigidities. The cumulative NNS distributions are plotted in 
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Fig. p. One can see that the NNS distribution corresponding to the anti-unitary invariant 
case, <p = 7r/2, is close to, but does not reach completely the GOE behavior. Similarly, 
for the <p = 0.47T case, the distribution is very close, but does not reach the GUE one. 
Nevertheless, the symmetry breaking is apparent and the numerical spectrum at ip = OAtt 
shows much more level repulsion than the GOE case, which is a clear-cut signature of the 
breaking of any anti-unitary symmetry. 

To measure departures from the entirely ergodic behavior, we can also use the Berry- 
Robnik distribution. The Berry-Robnik model for fully broken antiunitary symmetry con- 
sists of the superposition of two independent Poisson and GUE spectra. The results of the 
fits are collected in Table [IT]. We have not applied the Brody distribution to the broken 
anti-unitary invariance case, as the distribution is defined only for Brody parameter less 
than unity, and does not make any sense in the unitary case. The Izrailev distribution does 
not suffer this severe problem, and contain all limiting cases ((3 = for Poisson statistics, 
(3 = 1 for GOE and (3 = 2 for GUE). We have thus fitted our results with this distribution 
too. 

In all cases, the Izrailev distribution works much better than the Berry-Robnik ansatz, 
presumably because chaotic motion occupies most of the phase space and regular regions 
are very tiny as seen from Fig. [I] and from the values of q obtained. 

The values of the fitted Berry-Robnik parameter are consistent with the character of 
the corresponding classical motion. An independent information about the relative measure 
of a chaotic part of phase space comes from the fit of the A3 statistics, Fig. || [for broken 
anti-unitary invariance case A^ ^ is substituted by Af UE in Eq. (|3.2j )l, which turns out to 
agree very nicely with the values of the Berry-Robnik parameter. 

By gradually decreasing ip, we may observe the partial symmetry breaking by studying 
the variation of fitted parameters with ip. Such a transition has been analytically studied for 
gaussian random matrices [56| , |57| where the two-point correlation function was analytically 
found for an appropriate ensemble which interpolated between the GOE and the GUE. A 
further link with the dynamics of fully chaotic systems with partially broken antiunitary 



symmetry was also established ||58|| . These developments cannot be used here since the 
dynamics in our case is not fully chaotic (as seen from SOS plots and the non-integer level 
repulsion (3 parameters for extreme cases of preserved and broken antiunitary symmetry, see 
Table II). 

On the other hand, the Izrailev distribution is quite suitable since it should be a reason- 
able approximation both for a partial symmetry breaking and a mixed dynamics (it would 
be possible to construct an analog of Berry-Robnik distribution for such a case but it would 
be of little practical importance). 

Fig. [1^ summarizes the changes of the fitted Izrailev parameter (3 (small s repulsion) 
with (p. For ip = 7r/2, it is minimal and equal to 0.85 (see also Table [□]) for no around 55. 
With departure from the value <p = n/2 where the anti- unitary symmetry exists, it rapidly 
increases (filled circles) up to the maximal value 1.73 for ip = OAtc, 0.6tt. For still lower 
(higher) values of tp, the trend is reversed and (3 starts to decrease. This, at first glance, is 
a surprising effect (since the symmetry should not yet start to restore). However, it is most 
probably due the fact that the classical motion becomes more regular as <p is far from tt/2. 
Observe in Fig. [7] that the SOS around <p = is much more regular than for ip = tt/2. 

In order to test this hypothesis, we have used another n value. Indeed, if the dip of (3 
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near p = n/2 was due to classical reasons, it should not depend on % ,i.e. n . Fig. [10| shows 
also the fitted (3 parameters for n around 100. Generally, the (3 values obtained are slightly 
larger than for lower n . Clearly, p starts to decrease again around the value p> = OAir. On 
the other hand, the dip of (3 when tp goes to tt/2 (where an anti-unitary symmetry exists) 
is faster than for lower no- This is in a full qualitative agreement with RMT 0,[J: for the 
GOE — > GUE transition the parameter controlling the transition is proportional to iV -1 / 2 
where N is the matrix dimension. The size of our matrices scales as n 2 ,, so the parameter 
controlling the breaking (i.e. a deviation from 7r/2 value) should scale like 1 oc % e s- Such 
a behaviour is roughly observed in Fig. 10. 



Fig. 11 shows the "maximal" repulsion case obtained, i.e. data for n around 100 and 



<p = OAir. The numerical data are presented in the form of the histogram of spacings 
and compared with the Izrailev distribution (the fit has been performed, as usual, for the 
cumulative distribution; the resulting (3 = 1.83 value has been used to plot the Izrailev 
distribution). The dash-dotted and dashed curves correspond to Wigner GOE and GUE 
distributions, respectively. The fact that we observe level repulsion much stronger than the 
GOE behaviour is a signature of anti-unitary symmetry breaking. 

As mentioned before, our first results on the manifestations of symmetry breaking in our 
system have been obtained for the microwave frequency being twice the Kepler frequency, i.e. 
for the 2:1 resonant driving ||31|| . This choice was motivated by SOS plots in the absence of 
the electric field - see Fig. [I] - showing smaller regular islands for higher microwave frequency. 
However, as we have seen above, the presence of the electric field makes the secular motion 
in the principal resonance island chaotic enough and in fact we get stronger repulsion for the 



1 resonance than for the 2 : 1 situation reported before |31|. For completeness, we show 



in Table III the fitted parameters obtained from the numerical results reported in [31] using 



either the spacing distribution or the spectral rigidity A 3 . As can be seen, the conclusions 
we obtain from these results completely confirm the analysis of the 1:1 resonance. 

As far as we know, the studied system constitutes the first experimentally realizable 
example of a quantum system with broken anti-unitary symmetry. We have considered small, 
but finite, field amplitudes to stay well within the applicability range of the perturbation 
theory. Nevertheless, this is experimentally feasible: for n ~ 55 and F ~ 5 x 10~ 4 , i.e. 
about 0.3 V/cm, the mean level spacing is of the order of MHz. 

For stronger fields, our classical studies also suggest a similar behavior. It could even 
be that the breaking of the secular approximation makes the system more chaotic and that 
the statistical properties will be closer to GUE. However, we are not able to show quantum 
numerical results as they require full quantum numerical treatment which is difficult with 
the present computer resources and must be left for a future work. 



V. CONCLUSIONS 

We have considered an hydrogen atom perturbed by a resonant elliptically polarized 
microwave field with or without an additional static electric field of different orientation, in 
the limit of small field amplitudes. Classically, such fields may produce chaotic dynamics in 
the secular motion of the electronic elliptical trajectory. In quantum mechanical language, 
states coming from a given hydrogenic manifold may be mixed significantly only with each 
other. Such a situation has been interpreted as a signature of an intramanifold chaos. 
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For the pure microwave problem, we have studied two different resonant driving cases, 
i.e. the 1:1 and 2:1 resonances between the microwave field and the unperturbed Kepler mo- 
tion. Quantizing the fast orbital electronic motion, one can derive an effective Hamiltonian 
describing the slow secular precession of the electronic elliptical trajectory. For a generic 
elliptical polarization, the effective Hamiltonian has two degrees of freedom and turns out 
to be non-integrable. By means of Poincare surfaces of sections, we have found a range of 
the secular energy where the phase space reveals mixed character with a significant amount 
of chaotic layer. Switching to a quantum perturbation calculation, we have shown that the 
statistical properties of the corresponding quasienergy levels reveal an intermediate behavior 
between Poisson and GOE character. For the 2:1 resonance case, the classical phase space is 
significantly more irregular than for the 1:1 case and, consequently, the spectral properties 
are closer to the GOE behavior. 

The application of an additional static electric field to the system has allowed us to 
enhance chaos in the secular motion. Moreover, the static electric field, for a generic orien- 
tation, breaks any anti-unitary symmetry of the system which has a dramatic effect on sta- 
tistical properties of quasienergy levels. This is the first, to our knowledge, experimentally 
realizable quantum system, with corresponding chaotic classical behavior, which exhibits 
breaking of any generalized time-reversal symmetry. 

We have studied the principal 1:1 resonance for two slightly different static field ori- 
entations: the first one corresponds to the anti-unitary invariance case, the other one to 
breaking such a symmetry. The classical phase space structures, in both similar 
with predominately chaotic behavior. However, the statistical properties of the quantum 
spectrum change from a near-GOE to a near-GUE behavior when one switches from the 
preserved to broken anti-unitary symmetry case. In the intermediate situations, we could 
observe a partial symmetry breaking effect due to the finite size of matrices involved in the 
problem. 

We have studied the limit of small field amplitudes as it allows us to employ quantum 
perturbation theory. For higher amplitudes, the classical effective Hamiltonian is known, but 
finding the quasienergy spectrum requires full quantum numerical calculations. For a pure 
microwave perturbation, however, we may predict that statistical properties of the quantum 
spectrum, for stronger field, should be the same as for weak field limit. This is because 
of the specific form of the effective Hamiltonian, Eq. ( |2.17| ), which depends on dynamical 
variables only through the T m . 
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VII. APPENDIX A 



The amplitude and the phase of the atomic dipole at the microwave frequency, 
T m (L, *, M, a) and S(L, *, M, $; a) appearing in Eq. ( |2.14| ) may be expressed as ||| 



{( 



1 + a 



2 

I -a 2 



V 2 + U 2 



+ 



a 



V 2 + U 2 



(V m V- m - U m U- m ) cos 2$ - (VmU-m + U m V- m ) sin2$] } 



1/2 



(7.1) 



and 



tan 5 



(1 - a) ( V_ m sin $ + £/_ m cos$) - (1 + a)V r m sin<J) + L/" m cos$) 
(1 - a)(l / _ m cos$ - £/_ m sin$) + (1 + a)KjCOS$ - f/ m sin$) ' 



(7.2) 



where and ?7 m are Fourier expansion terms of the original hamiltonian ( |2.1|) in action- 
angle variables. Explicitly, they read 



3e 

K,(J,L,*) = -yJ 2 cos*, 
U (J,L,M,V) = ~^~ J 2 sm*, 



and for m ^ 



(7.3) 



K»(J,L,M,*) = — [J m (me) + M ^ g2 J m {me)\ cos *, 



tU J, L, M, ¥) = -tf J» + ^— 

777 C 



■JmW] sin*. 



(7.4) 



In the above formulae e = y 1 — L 2 / J 2 is an eccentricity of the electronic ellipse while 
J m {x) and j' m {x) denote the ordinary Bessel function and its derivative, respectively. 



VIII. APPENDIX B 



We give here explicit expressions for various level spacing distributions which have been 
used in the present paper. 

The Poisson distribution, corresponding to a system with classically integrable dynamics 
@, reads 

P(s) = exp (s) . (8.1) 

For an ergodic classical behavior, the quantum spectrum of a generic system is conjectured 
to have a nearest neighbor spacing (NNS) distribution (for the unfolded spectrum) similar to 
that of the random matrices of the same universality class. The resulting NNS distributions 
are quite complicated (see e.g. 0,§]). However, a good approximation is given by the so 
called Wigner surmise, obtained for matrices of rank 2. These are: 
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TV / 7T o 

P(s) = -sexp (^--s 
for an anti-unitary invariant (GOE) system and 



32 



P(s) = —s exp 



7T^ 



7T 



(8.2) 



(8-3) 



for broken anti-unitary invariance (GUE). 

The phenomenological Brody distribution |52 which interpolates between Poisson and 
GOE distributions reads 



P( s ) = C(P + l)/exp (-Cs p+1 



with 



C 



'(3 + 2' 



•4) 



(8.5) 



(3 = (resp. (3 = 1) corresponds to the extreme case of the Poisson (resp. GOE) statistics. 
Another attempt towards a simple distribution interpolating between different ensembles, 



is due to Izrailev p3| , |59l and reads 

P(s) = As^l + sBpy^exp 

where 



n 2 Ps 2 7T 
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(8.6) 



/(/?) = 2/3(1 >^/ 2 ) _ 0.3GS7I. 



and A, B are constants that ensure 



and 







P(s)ds 



sP(s)ds 



•7) 



(8.8) 



(8.9) 



It is claimed to work reasonably well for all possible intermediate situations. 

While there exist several other propositions in the literature, we list only the so called 



Berry-Robnik statistics ||51|| . It may be derived assuming an independent superposition of 
Poisson spectrum and spectra corresponding to random matrix predictions. If one deals 
with Poisson and only one GOE spectrum, the Berry-Robnik distribution reads 



P(s) 



71 



2q{l-q) + -q 3 s 



exp 



7T 



+ (1 — q) exp[(g — l)s]erfc —z~qs 



5.10) 
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where < q < 1 is the relative weight of the GOE spectrum. Classically, q corresponds to the 
relative volume of the chaotic part of phase space. The similar Berry-Robnik distribution 
for a superposition of a Poisson spectrum and one GUE spectrum is 



P(s) 



?|gV + V(l-?)« + (l-9) s 



7T 



7T 



exp I q 2 s 2 



+ 



2q(l — q) — (1 — q) 2 qs erfc ( —j=qs 1 > exp [(q — l)s] . 



Ul) 



Let us present also the expression for the U(W) function which has been employed in 
a fine-scale representation of the deviation of the numerical level spacing distribution from 



the best fitting theoretical distribution. The following function [61 



U(W) 



arccos 



Vi - w, 



[8.12) 



where W is the value of the cumulative level spacing distribution / P(s)ds, ensures that, 
over the whole range of W, i.e. from to 1, the standard deviation of numerical data is 
uniform and equal to SU = l/(n\fN), where N is the total number of spacings. 

For completeness, let us finally define the spectral rigidity As(L) as an average over the 
spectral range used in analysis (i.e. over x ) of 



pxo+L 

A 3 0o, L) = L _1 min A ,B / dx(N(x) - Ax- B)\ 
where N(x) is the integrated level density (a staircase function). 



(8.13) 
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TABLES 



TABLE I. Fitted parameters for different distributions as denned in the Appendix. The data 
collected correspond to the hydrogen atom driven by an elliptically polarized microwave field and 
are shown in Fig. ||| and Fig. ||. For each set of levels, we fit the nearest neighbor spacing distribution 
or the spectral rigidity A3 and show here the value of the free parameter for the best fit. 

1:1 resonance 2:1 resonance 

n = 50 - 59 n = 97 - 102 n = 50 - 59 n = 97 - 102 

q (Berry-Robnik): 056 0A7 0/77 073 

q (A 3 ): 0.55 0.43 0.76 0.73 

(3 (Izrailev): 0.17 0.12 0.42 0.37 

P (Brody): 0.23 0.16 0.49 0.44 

Time-reversal invariance Yes Yes Yes Yes 



TABLE II. Fitted parameters for different distributions as defined in the Appendix. The data 
partially shown in Fig. || correspond to the hydrogen atom driven by an elliptically polarized 
microwave field and exposed to an additional static electric field. For each set of levels, we fit the 
nearest neighbor spacing distribution or the spectral rigidity A3 and show here the value of the 
free parameter for the best fit. 
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TABLE III. Fitted parameters for different distributions as defined in the Appendix. The 
data correspond to the 2:1 resonance of the hydrogen atom driven by an elliptically polarized 
microwave field and exposed to an additional static electric field. For each set of levels, we fit the 
nearest neighbor spacing distribution or the spectral rigidity A3 and show here the value of the 
free parameter for the best fit. The data with scaled energy 7i G (8.5, 9.5) and for T = 5000 are 
analyzed. 
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FIGURES 



FIG. 1. Poincare surfaces of sections (at $ = 0) of the classical secular motion, Eq. ( |3.1|) , 
for the hydrogen atom perturbed by a weak resonant microwave field with elliptical polarization 
(polarization parameter a = 0.4). The coordinates used for the plot are the scaled total angular 
momentum Lq = L/uq and its canonically conjugate angle Left column - the 1:1 resonance case 
(microwave frequency equal to the unperturbed Kepler frequency of the electron on its elliptical 
trajectory) for the secular energies (going from bottom to top), H eS /F$ =0.02, 0.03, 0.04, 0.07. 
Right column - the 2:1 resonance case (microwave frequency is twice the Kepler frequency), for 
H c{i /F3 =0.0006, 0.0013, 0.0018, 0.0021 again from bottom to top. Note that, for the parameters 
chosen, not the whole (Lq,^) space is accessible. 

FIG. 2. Energy levels of the no = 20 hydrogenic manifold resonantly driven by an elliptically 
polarized microwave field vs. the degree of polarization, a, for the 1:1 resonance case [panel (a)] 
and 2:1 resonance case [panel (b)]. We plot the shifts (divided by Fq) of the energy levels with 
respect to the unperturbed energy of the atom. The upper states evolve smoothly with almost 
exact level crossings, in agreement with the classical mostly regular dynamics. In contrast, most 
of the lowest states evolve irregularly and sometimes display large avoided crossings, a signature 
of a classically chaotic behaviour. 

FIG. 3. A part of the spectrum of no = 50 hydrogenic manifold resonantly driven by an 
elliptically polarized microwave field vs. the degree of polarization, a, for the 2:1 resonance case in 
the region of mixed regular and chaotic motion. We plot the shifts (divided by Fq) of the energy 
levels with respect to the unperturbed energy of the atom. Most of the states display large avoided 
crossings, a signature of a classically chaotic behaviour. 

FIG. 4. Cumulative level spacing distribution, W(s), for the pure microwave perturba- 
tion (degree of the polarization a = 0.4) for the 1:1 resonance case. Levels in the range of 
T^-cs/Fq = 0.02 — 0.045 are analyzed. In panel (a), the solid line represents numerical data for 
no = 50 — 59 (there are about 10000 spacings); the dashed and dash-dotted lines represent the 
Poisson and GOE distributions respectively. Panel (b) shows a fine-scale representation of the 
deviation of the numerical level spacing distribution from the best Izrailev distribution in terms 
of the U(W) — U (Wi zra _n ev ) vs. W; the transformation U{W) = arccos y/1 — W is used in order 
to have uniform statistical error over the plot - compare |6l[] and the Appendix. The upper and 
lower noisy curves represent one standard deviation from the calculated data which lie in the mid- 
dle of the band. The long-dashed curve represents the best Berry-Robnik distribution while the 
dotted line the best Brody distribution. Panel (c) , the solid line represents the numerical data 
for no = 97 — 102 (there are about 20000 spacings), the dashed and dash-dotted lines represent 
the Poisson and GOE distributions as in (a). Panel (d) shows the corresponding deviations of 
the numerical level spacing distribution from the best fits - the notation is the same as in panel 
(b). Note that the best fitted Berry-Robnik or Izrailev distributions are altogether in excellent 
agreement with the numerical results. At the scale of (a) and (c), they are not distinguishable 
from the data. 
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FIG. 5. The same as in Fig. ||, but for the 2:1 resonance case for the range of 
7i cS /F$ = 0.00035 - 0.0018. 



FIG. 6. Spectral rigidity A3 for the pure microwave perturbation (a = 0.4) compared with 
the random ensemble predictions. Numerical data (solid lines) for no = 50 — 59 are shown in panel 
(a) for the 1:1 resonance case and in panel (b) for the 2:1 resonance. The data for no = 97 — 102 
are presented in panels (c) and (d) for the 1:1 and 2:1 resonances respectively. The dotted lines 
are the fits of Eq. (3.2) while the dashed lines indicate the Poisson (upper straight lines) and GOE 
(lower curves) predictions. Note the saturation at large L, in agreement with the semiclassical 
prediction. 



FIG. 7. Poincare surfaces of section (at $ = 0) of the classical secular motion, Eq. ([4.2j), for the 
hydrogen atom in a static electric field and driven by an elliptically polarized (a = 0.4) microwave 
field resonant with the Kepler frequency (1:1 resonance). Left and right columns correspond to a 
static electric field with scaled amplitude T = 10 and different orientations 99 = 0, 8 = ir/4 and 
ip = vr/2, 9 = 7r/4 respectively. The secular energies (going from bottom to top, the same for left 
and right panels) are: 7i = —0.2, 0.4, 0.8, 1.4. Note that, for the parameters chosen, not the 
whole (Lq,^) space is accessible. 



FIG. 8. Cumulative level spacing distribution, W(s), for the the hydrogen atom in an ellipti- 
cally polarized microwave field (ellipticity parameter a = 0.4) combined with a static electric field, 
for the 1:1 resonance case. Levels in the range of Ti = — 1.4 and for no = 50 — 59 are analyzed 
(about 20 000 spacings). The value of the reduced microwave strength is chosen as T = 10, while 
the angle is 6 = ir/4. Panel (a) shows the data (solid line) for <p = 0.4-7T, i.e. for broken generalized 
time reversal invariance. The fitted Izrailev distribution cannot be distinguished from numerical 
data. Dashed and dash-dotted lines correspond to GUE and GOE predictions respectively (thus 
the data are closer to GUE than GOE). Panel (b) shows the difference between the numerical 
result and the best Izrailev distribution (horizontal line). The upper and lower noisy curves yield 
one standard deviation from the numerical data that lie in the middle of the band. The Izrailev 
distribution stays well within the one standard deviation practically everywhere. The dashed line 
corresponds to a Berry-Robnik distribution (for the mixture of Poisson and GUE spectra) which 
is clearly inferior to the Izrailev distribution. Panel (c) and (d) correspond to <p = it/ 2 where 
the anti-unitary symmetry is restored. In (c), dashed line corresponds now to Poisson, the data 
trace the distribution close to GOE one. In panel (d) the dashed-dotted line corresponds to Brody 
distribution. Clearly Brody and Izrailev distributions practically coincide and work much better 
than the Berry-Robnik distribution. 
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FIG. 9. Spectral rigidity, A3 compared with the random ensemble predictions for the hydrogen 
atom resonantly driven by an elliptically polarized microwave field in the presence of a static 
electric field of different orientation. The numerical data used are the same as in figure § Panel 
(a) corresponds to broken anti-unitary symmetry {<p = OAtt), panel (b) to ip = tt/2 (generalized 
time-reversal invariant case). The dotted lines denote the fit of Eq. ( |3.2[) for (b) and a similar 
expression for (a) with GUE in place of GOE. The dashed lines indicate the Poisson (upper straight 
lines), GOE and GUE predictions. Observe that, for broken anti-unitary symmetry, the data trace 
between the GOE and GUE curves. The typical saturation effects appear for large L. 

FIG. 10. The gradual symmetry breaking with the orientation of the electric field, i.e. ip. Filled 
circles correspond to the fitted Izrailev repulsion parameters (3 for data around no = 55, open circles 
for data around no = 100. Close to ip = 0.5-7T - where an anti-unitary symmetry exists - there is 
sudden drop in (3 by roughly one unit, as predicted by Random Matrix Theory. For smaller (larger 
- the figure is symmetric around vr/2) values of ip, the classical motion is less chaotic with large 
regular islands in the phase space - this explains the decrease of (3 for ip < 0.4-7r(> 0.67r). For a 
further discussion, see text. 

FIG. 11. Nearest neighbor spacing distribution for the hydrogen atom resonantly driven by 
an elliptically polarized microwave field in the presence of a static field oriented in the direction 
(p = OAtt, 9 = it/ 4), such that the generalized time-reversal symmetry is broken. The numer- 
ical data correspond to no = 97 — 102 (40 000 spacings), the solid line gives the fitted Izrailev 
distribution. The dash-dotted and dashed lines correspond to GOE and GUE Wigner surmises, 
respectively. The numerical results clearly show a stronger level repulsion than in the GOE case, 
a signature of the breaking of any anti-unitary symmetry. 
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